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ABSTRACT 

This paper presents a re- analysis of Lauer & Postman's (1994) finding that the Abell 
cluster inertial frame (ACIF), defined by the 119 Abell clusters within 15,000 km s _1 , 
is moving at almost 700 km s _1 with respect to the cosmic microwave background. 
Such a motion is inconsistent with most cosmological models at a confidence level of 
95% or higher. We examine the use of the relation between the metric luminosity of 
brightest cluster galaxies and the slope of their luminosity profiles as an estimator of 
distances and peculiar velocities. We obtain an exact expression for a cluster's peculiar 
velocity in terms of the residual magnitude about this relation and compare this to the 
approximation used by Lauer & Postman. We critically examine the method used by 
Lauer & Postman to recover the Local Group motion from the scatter in this relation, 
and develop improved procedures including a maximum likelihood method that 
provides a direct estimate of the uncertainty in the derived motion. Simulations show 
this method yields an unbiased estimate for the Local Group motion with significantly 
smaller uncertainties than LP's original method. We re- analyse Lauer & Postman's 
data to obtain an improved estimate for the motion of the Local Group. We find that 
the Local Group is moving relative to the ACIF at 626 ± 242 km s" 1 towards /=216°, 
b= — 28° (±20°). This implies that the ACIF is itself moving relative to the cosmic 
microwave background at 764 ± 160 km s -1 towards /=341°, 6=49° (±20°). This 
motion is consistent with that derived by LP but has a 10% larger amplitude and 20% 
smaller uncertainties, making it even harder to reconcile with cosmological models. 

Subject headings: galaxies: clustering — galaxies: distances and redshifts — galaxies: 
elliptical and lenticular, cD — galaxies: kinematics and dynamics — large-scale 
structure of Universe 
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1. INTRODUCTION 

Measuring the bulk motion of large volumes of the Universe has proved a powerful tool for 
testing cosmological models, since the amplitude of such motions is directly related to the power 
spectrum of density fluctuations (Peebles 1993). With the development of more precise distance 
estimators in recent years the size of the volumes over which peculiar velocities can be measured 
has steadily increased (Burstein 1990, Scaramella et al. 1994, Dekel 1994). Now Lauer & Postman 
(1994; hereafter LP) have measured the motion of the Local Group with respect to (w.r.t.) the 
inertial frame defined by the mean motion of the 119 Abell clusters closer than c,?=15,000 km s _1 . 
They find that the Local Group is moving at almost 600 km s _1 w.r.t. this Abell cluster inertial 
frame (ACIF) and in a direction that is 75° away from the direction in which the Local Group is 
moving w.r.t. the cosmic microwave background (CMB) frame (assuming that the CMB dipole is 
indeed due to relative motion). This implies that the ACIF is itself moving w.r.t. the CMB frame 
at almost 700 km s _1 . 

That so large a volume of the Universe should be moving at such a high velocity w.r.t. the 
cosmic rest frame defined by the CMB is an intuitively surprising result. It implies that the CMB 
dipole must be generated by mass fluctuations at distances beyond those of the clusters in the 
ACIF, i.e. greater than 150 h _1 Mpc. As recent calculations and simulations show (Feldman & 
Watkins 1994, Strauss et al. 1994), the power spectra of all the usual cosmological models predict 
that such a large motion for the ACIF has a probability of less than 5%. 

Testing the correctness of this very important result is clearly essential. Several programs to 
do this are in progress: Lauer and Postman, with M.A.Strauss and J.Huchra, have begun work 
extending the cluster sample out to 24,000 km s _1 , increasing the sample size as well as the scale 
on which the bulk motion can be measured; the EFAR project (Colless et al. 1993) will obtain 
bulk motions in two superclusters using the D n —a distance estimator for ellipticals; similar studies 
are also being carried out using the Tully-Fisher relation for spirals. However it is also necessary 
to confirm the methodology used to derive the bulk motions from such observations. This paper 
examines some of the methodological issues using LP's original dataset. 

LP use the M-a relation between the metric luminosity and the luminosity profile slope of 
brightest cluster galaxies (BCGs) as the estimator of distances and peculiar velocities for clusters 
in the ACIF. We show that LP's expression for a BCG's peculiar velocity in terms of its residual 
with respect to the M-a relation is an approximation, and derive the exact relation. We also 
examine the physical basis for LP's procedure of minimising the rms BCG peculiar velocity as 
a means of inferring the Local Group motion in the case where the magnitude residuals are 
dominated by the intrinsic errors. We conclude that in these circumstances a straightforward 
minimisation of the magnitude residuals themselves is to be preferred. Finally, we re-analyse 
LP's original dataset using three different (though closely-related) methods to recover the Local 
Group motion: a revised version of their scheme employing minimisation of the rms BCG peculiar 
velocity, direct minimisation of the magnitude residuals, and a maximum likelihood scheme that 
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allows estimation of confidence intervals without requiring Monte Carlo simulations of the data. 
We obtain a Local Group motion that is consistent with that derived by LP, but our improved 
methodology yields smaller uncertainties that increase the significance of this result. 



2. THE DISTANCE ESTIMATOR 



In this section we re-derive the distance (and peculiar velocity) estimator based on the use of 
the relation between metric magnitude and structure parameter a for brightest cluster galaxies 
(BCGs). 

We consider our distance and motion w.r.t. a cluster in the ACIF. The ACIF is defined by the 
mean motion of the clusters it comprises. The redshift z that we measure to one specific cluster 
may be decomposed as the product of three redshift components: 

(l + z) = (l + z ff )(l + z c )(l + z p ). (1) 

Here zjj is the 'Hubble redshift' that we would measure if neither we nor the cluster had any 
peculiar velocity w.r.t. the ACIF (i.e. zjj is the redshift corresponding to the cluster's true Hubble 
distance), z c corresponds to the peculiar velocity of the cluster w.r.t. the ACIF, and z, p corresponds 
to the peculiar motion of the Local Group w.r.t. the ACIF. Note that for a particular cluster with 
direction given by the unit vector r, z p is related to the Local Group peculiar velocity v p by 



where the minus sign results from the fact that if the Local Group motion is in the direction of 
the cluster then the cluster's redshift w.r.t. the Local Group is reduced. 

We want to use the absolute luminosity of the BCG inside some specified metric radius, 
L(ro), as a standard candle distance estimator for the cluster. The angle, (f>o, that r$ subtends for 
the observer is given by the angular diameter distance of the cluster, 

d A = r /<f> , (3) 

while the apparent luminosity within this angle, l((f>o), is given by the luminosity distance, 

These distances are related to the Hubble redshift of the cluster by 

(1a{zh) = (1 + z H )~ 2 d L {z H ) (5) 

and 



d L (zH) = -TT- 
n 



I + z H + (1 + 2q z H 



,1/2 



1 + qoz H + (1 + 2q z H y/ 2 



(6) 
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We wish to measure L(ro), which means measuring l((f>o). However we require some estimate 
for zjj in order to determine (f>Q. Since our peculiar velocities are defined w.r.t. the ACIF, which is 
in turn defined by the mean motion of the clusters it comprises, the mean value of z c is zero and we 
may reasonably assume that z c <C zjj- Thus for a given value of z, p we can use the approximation 

(l + z ff )«(l + z)/(l + z p ) = (l + z')- (7) 
Using z' instead of zjj corresponds to measuring the luminosity, l((f>), inside the angle 

<f> = r /d A (z') = r/d A (z H ). (8) 

Thus we measure the luminosity inside a metric radius r rather than tq. 

If around r$ the BCG's luminosity profile can be approximated by a power law of slope a, 
then the ratio of the luminosity we actually measure to the one we seek is 

L(r) ( r 



L(r ) \r . 

( d A (z H ] 



V d A (z') 

( d L {zu)\ a ( 1 + z' x 2a 



V d L (z') J \1 + z H 
f d L (z H )\ a n i ^ 

V d L (z 



'^) a (l + z c f» (9) 



Now in practice we measure an apparent magnitude m((f>), which we have to convert to an 
absolute magnitude M(r). Absolute and apparent magnitudes are related by 

M{r) = m{4>) - 5 \og{d L {z H )) - K{z) - A - 25 (10) 

where K(z) is the K-correction for a galaxy at measured redshift z, and A is the absorption. Not 
knowing zjj, of course, we have to use z' , so we obtain 

M'(r) = m(4>)- 5 \og(d L (z')) - K(z) - A - 25 (11) 
= M(r)-51og(d L (z')/dL(*ff)). (12) 



The difference between this observed absolute magnitude and the true metric magnitude is 

M'(r) - M(r ) = M(r) - M(r ) - 5 \og(d L (z') / d L (z H )) 

= -2.5 log( L(r)/L(r j) - 5 log(d L (z')/d L (z H )) 

= -2.ba\og(d L (z H ) I d L (z')) - 5alog(l + z c ) - 5 \og(d L (z') / d L (z H )) 

= -2.5(2 - a) log(d L (z')/<*L(zff )) - 5a log(l + z c ). (13) 
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If we have assumed the correct Local Group motion then the above equation and the definition 
of z' tell us that M'(r) — M{tq) is only non-zero (in the absence of errors) if z c is non-zero. Since 
by definition the mean value of z c is zero, we expect that M'(r) will be very close to M{tq) in the 
mean. 

In fact in order to take advantage of the known correlation between M{tq) and a and thus 
have a better standard candle, we use as our estimate of the true metric magnitude the mean 
relation M'(r,a) evaluated at a((f>), assuming that this is a good estimator of M(ro,a). Thus the 
relation we use to estimate distances and peculiar velocities is 

AM = M'{r) - M'(r,a) = -2.5(2 - a) log(d L (z')/d L (z H )) - 5alog(l + z c ). (14) 

For an assumed Local Group motion v p and observed values of z, AM and a, we recover zjj and 
z c by simultaneously solving this equation and equation 1. 

This estimate of distance and peculiar velocity differs from that used by LP in two ways: 
(i) they omit the term 5alog(l + z c ) in equation 14, originating in the distinction between and 
d,L (see equation 9); (ii) they use the approximation 

d L (z')l d L (z H ) k, cz'/czh ~ 1 + cz c /cz H ~ 1 + cz c /cz (15) 

in deriving their distance estimator. (It should be noted, however, that LP do use the correct 
cosmological expressions for the angular diameter distance and luminosity distance in computing 
metric apertures and luminosities.) Making these changes to the expression for AM gives 

AM « -2.5(2 - a) log(l + cz c /cz). (16) 

from which LP obtain their expression for the cluster's peculiar velocity: 

" l) • (17) 

Figure 1 shows the difference, as a function of AM, between the exact value of cz c (obtained by 
solving equations 1 and 14) and LP's approximate value (given by equation 17) for an illustrative 
case. We assume i?o=80kms _1 Mpc -1 (though this does not enter the comparison) and go=0.5, 
and assume the cluster is at about the mean redshift of the ACIF sample, z=0.03, and the BCG 
has a typical structure parameter, a=0.6. We plot the relations for cz p = ±600 km s _1 , roughly 
the amplitude of the Local Group motion w.r.t. the ACIF found by LP. The top panel of the 
figure shows the value of cz c inferred from a given value of AM using the exact expression 
(cz p =+600 km s _1 solid line; z p = — 600 km s _1 dashed line) and LP's approximation (thick 
solid line; the two values of cz p give identical relations in LP's approximation). As the middle 
panel shows, the difference Acz c =cz c (LP) — cz c (true) can be large (several hundred km s _1 ) for 
values of AM that are well within the range found by LP (solid line shows the difference for 
cz p =+600 km s _1 , dashed line for cz p = — 600 km s _1 ). The bottom panel shows the fractional 
difference, which ranges up to 50% for observed values of AM. 
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Using LP's approximation to infer peculiar velocities for individual clusters can thus result in 
significant errors which are a function of the Local Group motion. However LP do not use their 
distance indicator for individual objects, but instead use it to derive the Local Group motion 
w.r.t. the ensemble of clusters making up the ACIF. By averaging over this all-sky sample the 
errors resulting from the use of the approximation cancel in first order, and the residual bias can 
be corrected for using Monte Carlo simulations. In the following section we examine the method 
used by LP to derive the Local Group motion and consider alternative approaches. 



3. ESTIMATING THE LOCAL GROUP PECULIAR VELOCITY 



The method LP use to derive the Local Group peculiar motion is to note that the rms peculiar 
velocity of the clusters in the ACIF will be minimised if it is computed using the true value of 
the Local Group motion. If an incorrect Local Group motion is assumed in computing the cluster 
peculiar velocities, this will add some spurious peculiar velocity Acz p to cz c , so that equation 17 
becomes 
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For the values of AM and a derived from this incorrect value of v p , LP find the value of Av p that 
minimises the rms cluster peculiar velocity v c , given by 



1/2 



(19) 



where the index i runs over all the ACIF clusters and the W{ are appropriate weights. They then 
update their estimate of v p with this Av p and recompute AM and a for each object accordingly. 
The whole process is then repeated until Av p becomes negligibly small and the true value of v p 
has been recovered. The weighting scheme preferred by LP has 



w t = (2 - a t ) 2 /zf, 



(20) 



which is intended to minimise the errors in v p by allowing for the fact that more compact galaxies 
(with a — ► 0) and more nearby galaxies (with Z{ — ► 0) are more sensitive to v,, 



The physical basis for this approach to recovering v p rests upon the assumption that the 
residuals AM about the mean M-a relation are attributable to the peculiar velocities of the 
clusters w.r.t. the ACIF. However the observed rms residual of <7m=0.24 mag, if interpreted in this 
way, would imply an implausible rms cluster peculiar velocity of about 1300 km s _1 . LP attempt 
to separate the intrinsic scatter in the M-a relation from that due to the cluster peculiar velocities 
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Fig. 1. — Comparison of the exact relation between cz c and AM with that obtained from LP's 
approximation. The comparison assumes the cluster is at z=0.03 and the BCG has a=0.6 and 
shows the relations for cz p =+600 km s _1 (solid lines) and —600 km s _1 (dashed lines): (a) cz c 
as a function of AM using the exact expression (thin lines) and LP's approximation (thick 
lines; the two values of cz p give identical relations in LP's approximation); (b) the difference 
Acz c =cz c (LF) — cz c (true); (c) the fractional difference Acz c /cz c . 
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by using the fact that the latter depends on redshift, and come to the unsurprising conclusion that 
the intrinsic scatter must greatly dominate the scatter due to cluster peculiar velocities. 

That being so, there is no physical reason for preferring to minimise v c rather than gm- As 
Figure 1 shows, AM and cz c are monotonically (and almost linearly) related. Hence minimising 
gm over v p should be very similar to minimising v c over v p , the only difference being that the 
relative weighting of clusters with different values of AM will differ to the extent that cz c is not 
linearly related to AM. It is worth noting that, for small AM (i.e. small cz c ), LP's expression for 
cz c (equation 17) reduces to 

fczAM\ , s 

cz c = -0.41ii 10 (— —) . (21) 

Applying the weighting of equation 19 in computing the rms cluster peculiar velocity gives 

1/2 



0.4c In 10 



E,-(2 - ai y/z 



(22) 



Since the observed redshifts Z{ are fixed and a 4 - changes very little with v p , minimising the 
rms cluster peculiar velocity w.r.t. v p using these weights is, for small AM, almost identical to 
minimising the unweighted rms magnitude residual. 

An alternative to minimising v c or gm is to use a maximum likelihood approach based on 
the observation that the residuals about the mean M-a relation have a Gaussian distribution 
(Postman & Lauer 1995). The probability of observing a given BCG's magnitude residual given 
an assumed v p is P 4 - = N(AM : 8M,gm) where AM is the observed residual about the mean 
M-alpha relation, SM is the predicted residual given by the RHS of equation 14, gm is the 
minimised scatter in the M-a relation and N(x : fj,, g) is the value at a; of a Gaussian of mean fj, 
and dispersion g. Computing SM requires z c , but we can make the simplifying approximation 
SM=0 (corresponding to z c =0) because the intrinsic scatter completely dominates the cluster 
peculiar velocities, so that SM <C gm and hence P 4 - ~ N(AM : 0, &m)- 

Given this probability for each galaxy in the sample we can recover the Local Group motion 
by maximising the log-likelihood for the entire sample, ln£ = E; m (-f«')? as a function of the 
assumed value of v p . The particular advantage of this method is that it allows us to estimate 
uncertainties using the property of the likelihood that 21n(£ max / £) is distributed approximately 
as xl, with the number of degrees of freedom v corresponding to the number of free parameters 
in the model. We can therefore use constant-^ 2 boundaries to obtain reliable confidence intervals 
(uncertainties) for the model parameters both individually and jointly without recourse to Monte 
Carlo simulations (see, e.g., Press et al. 1986, pp533-7). 

We can thus estimate the Local Group motion by minimising over v p any one of four different 
(though intimately connected) quantities: (a) the weighted rms cluster peculiar velocity, v c (LP), 
using LP's approximation; (b) the weighted rms cluster peculiar velocity, v c , using the exact 
expression; (c) the unweighted rms magnitude residual in the M-a relation, gm', and (d) the 
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negative log-likelihood of the magnitude residuals, — ln£. The details of the procedure are as 
follows: 

1. Adopt a Local Group peculiar velocity v p . For each cluster i in the ACIF compute 

z p i = —fi ■ v p /c where f 4 - is the unit vector in the direction of the cluster. Then compute 

z[ = (1 + Zi )/(1 + Zpi) - 1. 

2. Estimate the approximate angle subtended by the metric radius, ^ = ro/dA(z'j), and 
measure the apparent magnitude mi((f>) and cti((f>) = dlog L / d log (j)\ c f >t from the observed 
luminosity profile of the BCG. 

3. Compute M'(r) from equation 11. Fit the M'(r) versus a relation with a quadratic and 
obtain the residuals about this mean relation, AM; = M'(r) — M'(r, a 4 -). 

4. For method (a): Compute the notional rms peculiar velocity of the ACIF clusters using LP's 
approximation and weights (equations 17, 19 and 20). 

5. For method (b): Solve equations 1 and 14 simultaneously to derive zjji and z c i, and compute 
the notional rms peculiar velocity of the ACIF clusters from equations 19 and 20. 

6. For method (c): Compute the rms residual gm about the mean M-a relation. 

7. For method (d): Compute the log-likelihood ln£ = J2i ln(A(AM : 0,<7m)) 5 where Wm is the 
minimum rms magnitude residual. 

8. Repeat these steps, attempting to minimise cjm, v c or — ln£ as a function of v p . The value 
of v p that gives the minimum is the estimate of the Local Group peculiar velocity. 



One further difference between this procedure and that used by LP is that we update the a 4 - and 
AM; at every step in our minimisation rather than using a two-stage process of minimising v c over 
v p for a fixed set of a 4 -'s and AM 8 's and then updating these quantities for the new estimate of v p 
and repeating the process until it converges. Using the one-step procedure is conceptually neater, 
though it probably makes little difference to the final result. 

We now apply these various methods to LP's data for the BCGs and clusters in the ACIF in 
order to see what differences result in the inferred Local Group motion. 



4. RE- ANALYSING THE ACIF 



In re-analysing the ACIF with the procedures described above, we use the original data on 
the Abell cluster sample and the BCGs given in Lauer & Postman (1994; LP) and Postman & 
Lauer (1995). Since some processing of the observed luminosity profiles is required to derive the 
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values of a, M and AM for a given v p , the first step is to make a consistency check between our 
values and those obtained by LP. 

Table 3 of LP lists the a, M and AM of each BCG for three assumed Local Group motions: 

(i) case C, in which the ACIF is assumed to be at rest w.r.t. the CMB, so that the Local Group 
velocity w.r.t. the ACIF is the same as the Local Group velocity w.r.t. the CMB, i.e. v p =C; 

(ii) case L, in which the Local Group and ACIF are assumed to be at rest w.r.t. each other, i.e. 
v p =0; and (iii) case F, in which the Local Group has LP's derived velocity L w.r.t. the ACIF, i.e. 
v p =L. (Note that case L does not correspond to v p =L and that case F does not correspond to 
v p =F, the motion of the ACIF w.r.t. the CMB.) There is a minor ambiguity in recovering LP's 
results due to the fact that there are 4 clusters (407, 419, 1177 and 3574) for which the cluster 
redshift (computed using the value in Table 1 of Postman & Lauer (1995) transformed to the 
Local Group frame as specified in LP) is not identical to the value listed in Table 1 of LP. These 
differences are small (the largest is less than 100 km s _1 ) and we have simply adopted the values 
from Postman & Lauer (1995) since they give somewhat better agreement with the a, M and AM 
values listed in Table 3 of LP. 

Figure 2 shows the differences between our values of a, M and AM and LP's values for their 
three tabulated cases of v p . (Note that the values given by LP for case F are for their estimated 
Local Group velocity without the small correction for the bias introduced by the sample geometry 
which they apply to obtain their final L.) Since we are attempting to reproduce LP's procedure 
for deriving these quantities, this is purely a test of whether we are successful in recovering 
their numbers, and is not a check on whether the assumed Local Group motion is correct. The 
differences are plotted as a function of cos 9, the cosine of the angle between the direction of the 
cluster and the direction of the assumed Local Group motion. For all three cases any systematic 
differences are negligible and the rms differences between our values and LP's are approximately 
0.003 in a, 0.001 mag in M and 0.004 mag in AM, corresponding to an rms difference in cz c of 
~25 km s _1 . 

This comparison reassures us that we are recovering a, M and AM in a manner consistent 
with LP, and we can proceed to recover v p using the methods described in the previous section. 
Table 1 gives the value of v p obtained using the four methods described above and starting 
from three initial estimates: Vi n a=0, Vi n n=C and Vi n n = — C. We used the AMOEBA routine from 
Numerical Recipes to do the minimisation, and demanded that the quantity being minimised 
converge to within a fractional tolerance of 10~ 4 . The table lists both the initial velocity Vi n n and 
the recovered Local Group motion v p in Cartesian Galactic coordinates. For this value of v p , the 
table gives the rms residual about the quadratic fit to the M-a relation (ctm), the notional rms 
cluster peculiar velocity (v c ) and the log-likelihood of the magnitude residuals (ln£). Note that 
only one of these three quantities is actually optimised at a time; the others are simply the values 
at the optimum v p . The final column is the confidence level, P(C), at which we can reject the 
solution based on it's likelihood ratio w.r.t. the maximum likelihood solution — this is a measure of 
the consistency of the solutions obtained from the different methods. 
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Fig. 2. — The differences between our derived vaiues of a, M and AM and those tabuiated in 
Tabie 3 of LP for three specific Local Group velocities: v p =C, v p =0 and v p =L. The differences 
are plotted as a function of cos 9, the cosine of the angle between the direction of the cluster and 
the direction of the assumed Local Group motion. The rms scatter is given excluding the four 
clusters (indicated by crosses) with slightly discrepant redshift values (see text). 
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TABLE 1 

The Local Group Peculiar Velocity 
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-0.2302 


<1% 


( +7, -542, +302) 


(-490, 


-373, 


-342) 


0.242 


1270 


-0.2364 


<1% 


( -7, +542, -302) 


(-466, 


-364, 


-343) 


0.242 


1270 


-0.2570 


<1% 


(d) Maximising ln£: 
















( o, 0, 0) 


(-506, 
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-329) 
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-0.2288 




( +7, -542, +302) 


(-508, 


-373, 


-328) 


0.242 


1268 


-0.2290 




( -7, +542, -302) 


(-507, 


-371, 


-330) 


0.242 


1269 


-0.2288 





TABLE 2 
Comparison with Lauer & Postman 





v x 


Vy 


v z 


\ v p\ 


(/p, 6p) 




Ad 


P(C) 




( km s _1 ) 


( km s _1 ) 


( km s _1 ) 


( km s _1 ) 


(deg) 


(deg) 


(deg) 




v p 


-506 ± 188 


-367 + 213 


-329 ± 163 


626 ± 242 


(216,-28) 


20 






L 


-470 ± 250 


-399 ± 273 


-333 ± 198 


561 ± 284 


(220,-28) 


27 


4 


0.8% 


C 


+7+ 8 


-542 + 15 


+302 + 15 


621+ 10 


(271, +29) 


1 


78 


>99.99% 























98.5% 
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For each method the final value of v p is similar for all values of V{ n n\ other values of Vi n a, not 
shown in the table, were also tried and gave the same values of v p . As expected from the above 
discussion, the values of P(C) show that all four methods give consistent results for v p within the 
uncertainties, although there is a slight systematic difference between the methods minimising the 
cluster peculiar velocities and those minimising the magnitude residuals. 

We estimate uncertainties from the maximum likelihood method. The 68% confidence 
interval for each component of v p =(V x ,V y ,V z ) is the projection of the 3D likelihood contour 
2(ln C max — In £) = 1.0 onto that component's axis. We find that the Cartesian components of the 
best solution for v p are 

V x = -506 ± 188 km s" 1 

V y = -367 ±213 km s" 1 (23) 
V z = -329 ± 163 km s" 1 

which corresponds to a Local Group motion \v p \ = 626 ± 242 km s _1 towards 
(lp,b p ) = (216,-28) with a total angular error ag = 20°. Note that the radial 
amplitude of the Local Group motion has been corrected for error biasing according to 
\ Vp \ 2 = (V 2 + Vy 2 + V z 2 ) - (6V* + SV 2 + SV Z 2 ) = 706 2 - 328 2 . Figure 3 shows how the rms cluster 
peculiar velocity v c , the rms residual magnitude gm and the log-likelihood ln(£) vary with the 
assumed amplitude \v p \ of the Local Group motion, while Figure 4 shows how these quantities 
vary with the assumed direction (l p ,b p ) of the Local Group motion. The dipole signature of the 
Local Group motion is clearly seen in Figure 4; the extent to which the pole and anti-pole are not 
180° opposed is consistent with the error estimates. The M-a relation for our best estimate of v p 
is shown in Figure 5. The best-fit quadratic to this relation has an rms residual of 0.242 mag and 
coefficients cq = — 20.856, c\ =— 4.529 and C2=2.845. We confirm the finding of Postman & Lauer 
(1995) that there is no significant variation of the residuals with a or z, and that the distribution 
of residuals is well-represented by a Gaussian, as we assumed in constructing the likelihood. 

Table 2 compares our maximum likelihood solution for v p with the preferred solution of LP 
(v p =L), the case in which the ACIF is at rest w.r.t. the CMB (v p =C) and the case in which the 
Local Group is at rest w.r.t. the ACIF (v p =0). We list the components of these velocities (and 
their uncertainties) in both Cartesian and spherical Galactic coordinates. We also give the angle 
between our value of v p and the alternate values, A0, and the confidence level, P(C), at which we 
can reject the alternate values based on their likelihood ratio w.r.t. our preferred solution. 

The solution we obtain for v p is entirely consistent with that of LP — the offset is just 
L — v p = (+36, —32, —4). Differences in the procedures for deriving v p have not led to any 
significant differences in the final result. However it is worth pointing out that the agreement 
in \v p \ is slightly worse than might be expected from the very good agreement in (V x ,V y ,V z ). 
This is entirely a consequence of the fact that we obtain roughly 20% smaller uncertainties for 
the components of v p and so make a smaller correction for error biasing in computing \v p \. 
These smaller uncertainties mean that we rule out the alternatives v p =C and v p =0 with greater 
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Fig. 3. — (a) The rms cluster peculiar velocity, (b) the rms residual magnitude, and (c) the 
log-likelihood, all as functions of \v p \ for an assumed direction of the Local Group motion of 
(/ p ,6 p )=(216,-28). 
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(a) RMS cluster peculiar velocity with direction of Local Group motion for |v LG |=600 km/s 




b=-90 



(b) RMS magnitude residual with direction of Local Group motion for |v LG |=600 km/s 

b=+90 




b=-90 



(c) Log-likelihood with direction of Local Group motion for |v LG |=600 km/s 

b=+90 




b=-90 



Fig. 4. — Contour plots of (a) the rms cluster peculiar velocity, (b) the rms residual magnitude, 
and (c) the log-likelihood, as functions of (l p ,b p ) for an assumed amplitude of the Local Group 
motion of |u p |=600 km s _1 . The inferred direction of motion of the Local Group w.r.t. the ACIF 
is indicated as 0; the direction of motion of the Local Group w.r.t. the CMB is indicated as ©. 
The dots show the positions of the ACIF cluster sample. The zone of avoidance lies between the 
dashed lines at 6=±15°. 
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Fig. 5. — The M-a relation with the best-fitting quadratic shown as the solid line and the ±1(7 
error range shown by the dotted lines. 
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confidence than LP and provide tighter constraints on cosmological models seeking to explain 
the amplitude of the observed bulk motion of the ACIF. In the next section we use Monte Carlo 
simulations to confirm these error estimates and also to check for biases in the various methods 
for deriving v p . 

Our result for the motion of the ACIF w.r.t. the CMB, F=C-v p , is 

V x = +513 ± 188 km s" 1 

V y = -175 ± 214 km s" 1 (24) 
V z = +631 ± 164 km s" 1 

which is of course also consistent with LP's result, though with 20% smaller uncertainties in each 
velocity component. This corresponds to |.F| = 764 ± 160 in the direction (/, 6) = (341,49) with a 
total angular error gq = 20°. 

5. MONTE CARLO SIMULATIONS 

We have constructed Monte Carlo simulations in order to test and compare the various 
methods for recovering the Local Group motion. The inputs to these simulations are the Local 
Group motion v p and the M-a relation and its intrinsic scatter gm- From these and the observed 
directions, redshifts and a's of the BCGs we construct a simulated set of observed magnitudes. 
We then run our procedure to recover v p from these simulated observations. The assumptions 
that we are making are: (i) since we use fixed a's, that these change little with different assumed 
values for v p (as is in fact the case); (ii) that the distribution of the magnitude residuals really is 
Gaussian; and (iii) that the clusters themselves have negligible (or at least uncorrelated) peculiar 
motions. 

For each method of recovering v p we ran 5000 simulations using as inputs the best estimates 
of v p , gm and the M-a relation given in the previous section. Figure 6 shows the distributions 
of the components of the recovered v p and the distribution of the inferred gm for each method. 
We find that when we maximize ln(£) or minimize gm there is no significant bias in the recovered 
v p — the amplitude of the difference between the input v p and the mean v p recovered from the 
simulations is about 6 km s _1 , less than 1.5 times the joint standard error in the mean. When we 
minimize v c , however, this difference is (+20,-33,-5), the amplitude of which is almost 9 times 
the joint standard error. The small but significant bias for this method results from the fact that 
one is mis-interpreting the magnitude residuals as peculiar velocities. Since the bias correction we 
derive for our preferred maximum likelihood method is both small and statistically insignificant, 
no such correction has been applied to the results given in the previous section. 

It is interesting to note that the recovered value of the intrinsic scatter in the M-a relation, 
gm, is biased low by 0.006 mag. This reflects the fact that we are obliged to use the mean M-a 
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Fig. 6. — The components of the Local Group motion v p and the scatter in the M-a relation gm 
recovered from 5000 simulations by (a) maximising ln£, (b) minimising njf, and (c) minimising 
v c . The input values for the simulations were v p = ( — 506, — 367, — 329) and <7m=0.242 mag. The 
numbers in parentheses in each panel are the mean and standard deviation from the simulations; 
the smooth curve is the corresponding Gaussian. 
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relation calculated from the data rather than the unknown true relation, and hence tend to 
under-estimate the scatter. Thus our estimated scatter of 0.242 mag implies that the intrinsic 
scatter in the M-a relation is in fact 0.248 mag. 

As well as showing us that the maximum likelihood and minimum gm methods are effectively 
unbiased, the simulations also demonstrate the reliability of the maximum likelihood estimates 
for the uncertainties. The standard deviations in the velocity components recovered from 
the simulations, (±180, ±197, ±144), are in good agreement with the uncertainties estimated 
by maximum likelihood, (±188, ±213, ±163). The slightly lower Monte Carlo errors are the 
consequence of the fact that in constructing the simulations we used the observed scatter of 
0.242 mag rather than the true instrinsic scatter of 0.248 mag. 

6. CONCLUSIONS 

We have re-examined Lauer & Postman's (1994; LP) analysis of the Local Group motion 
w.r.t. the inertial frame defined by the mean motion of all the Abell clusters within 15,000 km s _1 
(the ACIF). Particular attention has been paid to the use of the M-a relation for brightest cluster 
galaxies (BCGs) as the distance estimator, and to the algorithm used to recover the peculiar 
motion of the Local Group w.r.t. the ACIF from these distances. 

We derive an exact expression for the clusters' peculiar motions in terms of their magnitude 
residuals AM about the M-a relation. This result significantly improves on the approximation 
used by LP. We also show that, with the weighting scheme preferred by LP, obtaining an estimate 
of the Local Group motion by minimising the rms cluster peculiar velocity v c is closely related to 
a more straightforward minimisation of the rms residual about the M-a relation, gm- We argue 
that when the intrinsic scatter in the M-a relation dominates the scatter due to the clusters' 
peculiar velocities (as is the case here) then there is no physical basis for preferring to minimise 
v c rather than gm- In fact, mis-interpreting the magnitude residuals as due to peculiar velocities 
rather than intrinsic scatter results in a small bias in the derived Local Group motion. Using 
the observation that the magnitude residuals have a Gaussian distribution with intrinsic scatter 
(7^=0.242 mag, we construct a maximum likelihood method for estimating the Local Group 
motion v p . This method has the advantage over the minimisation approaches that it provides 
direct estimates of the uncertainty in v p and the confidence levels at which alternate values for v p 
can be rejected, without recourse to Monte Carlo simulations. 

We apply these methods to LP's original data to obtain an improved estimate of the Local 
Group motion. We find that the Local Group is moving w.r.t. the ACIF at 626 ± 242 km s _1 
towards /=216°, 6=-28° (±20°). This implies that the ACIF is moving w.r.t. the CMB frame at 
764 ± 160 km s _1 towards /=341°, 6=49° (±20°). These results are consistent with those of LP 
but have 20% smaller uncertainties. Calculations by Feldman & Watkins (1994) and Strauss et al. 
(1994) suggest that LP's original result for the bulk motion of the ACIF is inconsistent with most 
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current cosmological models at a confidence level of 95% or higher. The slightly larger amplitude 
and significantly smaller uncertainties that we derive here for the motion of the ACIF w.r.t. the 
CMB further tighten the constraints on cosmological models. 

A number of stimulating discussions with Alister Graham are gratefully acknowledged. Marc 
Postman provided a careful reading of the manuscript and his comments resulted in several 
material improvements. 
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